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Abstract 

When considering perturbations in an open (Qq < 1) universe, cosmologists retain only 
sub-curvature modes (defined as eigenfunctions of the Laplacian whose eigenvalue is less 
than —1 in units of the curvature scale, in contrast with the super-curvature modes whose 
eigenvalue is between —1 and 0). Mathematicians have known for almost half a century that 
all modes must be included to generate the most general homogeneous Caussian random 
field, despite the fact that any square integrable function can be generated using only the 
sub-curvature modes. The former mathematical object, not the latter, is the relevant one 
for physical applications. The mathematics is here explained in a language accessible to 
physicists. Then it is pointed out that if the perturbations originate as a vacuum fluctuation 
of a scalar field there will be no super-curvature modes in nature. Finally the effect on the 
cmb of any super-curvature contribution is considered, which generalizes to $7o < 1 the 
analysis given by Grishchuk and Zeldovich in 1978. A formula is given, which is used to 
estimate the effect. In contrast with the case r^o = 1, the effect contributes to all multipoles, 
not just to the quadrupole. It is important to find out whether it has the same / dependence 
as the data, by evaluating the formula numerically. 

1 Introduction 

On grounds of simplicity, the present energy density ^Iq of the universe is generally assumed 
to be equal to unity (working as usual in units of the critical density) It is not however 
well determined by observation [||]. The density of baryonic matter can only be of order 0.1 
or there will be a conflict with the nucleosynthesis calculation, and although non-baryonic 
matter seems to be required by observation ||2| there is no guarantee that it will bring the 
total up to J7o = 1- Nor should one assume that a cosmological constant or other exotic 
contribution to the energy density will play this role. 

From a theoretical viewpoint the value = 1 is the most natural, because any other 
value of is time dependent. The preference for = 1 is sharpened if, as is widely believed, 
the hot big bang is preceded by an era of inflation. In that case Q has its present value at 
the epoch when the present Hubble scale leaves the horizon, and for a generic choice of the 
inflaton potential this indeed implies that is very close to 1 more or less independently of 

^Throughout this article = 1 will mean a value of Qq close to 1, and f2o < 1 will mean a value 
substantially less than 1, say less than 0.9. 



the initial value of $7. It is also easier for inflation to explain the homogeneity and isotropy 
of the observable universe if f^o = 1- On the other hand it is certainly not the case that 
r^o = 1 is an unambiguous prediction of inflation j^, . 

The literature on the Oq < 1 cosmology is small compared with the enormous output 
on the case Qq = 1, because the latter is simpler and observations that can distinguish the 
two are only now becoming available. This is especially true in regard to the subject of the 
present paper, which is the effect of spatial curvature on cosmological perturbations. The 
only data relevant to this subject are the lowest few multipoles of the cosmic microwave 
background (cmb) anisotropy, that were measured recently by the COBE satellite (5|, 0. 

This article is concerned both with the basic formalism that one should use in describing 
cosmological perturbations, and with the cmb multipoles. To describe its contents, let us be- 
gin by recalling the presently accepted framework within which cosmological perturbations 
are discussed. 

Cosmological perturbations are expanded in a series of eigenfunctions of the Laplacian 
for two separate reasons. One is that each mode (each term in the series) evolves indepen- 
dently with time, which makes it easier to evolve a given initial perturbation forward in 
time. The other is that by assigning a Gaussian probability distribution to the amplitude 
of each mode, one can generate a homogeneous Gaussian random field. Such a field consists 
of an ensemble of possible perturbations, and it is supposed that the perturbation seen in 
the observable universe is a typical member of the ensemble. The stochastic properties of 
a Gaussian random field are determined by its two point correlation function (/(l)/(2)), 
where / is the perturbation and the brackets denote the ensemble average, and the adjective 
'homogeneous' indicates that the correlation function depends only on the distance between 
the two points. 

The question arises which eigenfunctions to use, and in particular what range of eigen- 
values to include. If f^o = 1 space is flat and it is known that the Fourier expansion, which 
includes all negative eigenvalues, is the correct choice. It is complete in two distinct re- 
spects. First, it gives the most general square integrable function, so that initial conditions 
in a finite region of the universe can be evolved forward in time. Secondly, it gives the most 
general homogeneous Gaussian random field. Instead of the Fourier expansion one can use 
the entirely equivalent expansion in spherical polar coordinates. 

If < 1, the curvature of space defines a length scale. The spherical coordinate 
expansion can still be used, and it is known [^, that the modes which have real negative 
eigenvalue less than —1 in units of the curvature scale provide a complete orthonormal 
basis for square integrable functions. Presumably for this reason, only these modes have 
been retained by cosmologists. We will call them sub-curvature modes, because they vary 
significantly on a scale which is less than the curvature scale. The other modes, with 
eigenvalues between —1 and in units of the curvature scale, we will call super-curvature 
modes. 

It is certainly enough to retain only sub-curvature modes if all one wishes to do is 
to track the evolution of a given initial perturbation, since the region of interest is always 
going to be finite and any function defined in a finite region can be expanded in terms of the 
sub-curvature modes. (In fact, to describe the observations that we can make it is enough 
to specify initial conditions within our past light cone.) But this is not what one does in 
cosmology.0 Rather, one uses the mode expansion to generated a Gaussian perturbation, 
by assigning a Gaussian probability distribution to the amplitude of each mode. In this 
context the inclusion of only sub-curvature modes looks restrictive. For example, it leads to 
a correlation function which necessarily becomes small at distances much bigger than the 

^The only case where one is interested in evolving a given initial condition is when one obtains the well 
known relation between the matter density contrast and the galaxy peculiar velocity field, but one uses 
this prediction only on very small scales where the curvature cannot be significant. Even then, one is still 
interested in the stochastic properties as well. 



curvature scale (to be precise, it is less than r j sinhr times its value at r = 0, where r is 
the distance in curvature units). 

Faced with this situation, we queried the assumption that only sub-curvature modes 
should be included, and the results of our investigation are reported here. 

First we describe the mathematical situation, showing that indeed a more general Gaus- 
sian random field is generated by including also the super-curvature modes. As expected 
the correlation function can now be constant out to arbitrarily large distances. 

Then we go on to ask whether nature has chosen to use the super-curvature modes, 
focussing on the low multipoles of the cmb anisotropy which are the only relevant observa- 
tional data, and on the curvature perturbation which is thought to be responsible for these 
multipoles. If, as is usually supposed, this perturbation originates as a vacuum fluctuation 
of the inflaton field, there will be no super-curvature modes. On the other hand, like any 
other statement about the universe one expects this assumption to be at best approximately 
valid. Supposing that it fails badly on some very large scale, hut that the curvature per- 
turbation still corresponds to a typical realization of a homogeneous Gaussian random field, 
one is lead to ask if a failure of the assumption could be detected by observing the cmb 
anisotropy. We note that for Qq = 1 this question has already been discussed by Grishchuk 
and Zeldovich |10|, and we extend their discussion to the case < 1- 

After our investigation was complete, and the draft of this paper was almost complete, 
M. Sasaki suggested to one of us (DHL) that a mathematics paper written by Yaglom in 
1961 might be relevant. From this paper we learned that the need to include both 
sub- and super-curvature modes in the expansion of a homogeneous Gaussian random field 
in negatively curved space has been known to mathematicians since at least 1949 [||]. It 
would appear therefore that the assumption by cosmologists that only the sub-curvature 
modes are needed is a result of a complete failure of communication between the worlds 
of mathematics and science, which has persisted for many decades. We have retained the 
mathematics part of our paper because it gives the relevant results in the sort of language 
that is familiar to physicists, though it is strictly speaking redundant. 

Let us end this introduction by saying a bit more about the cosmology literature. Start- 
ing with the paper of Lifshitz in 1946 |jl^, there are many papers on the treatment of 
cosmological perturbations for the case Q < 1. However, most of them deal with the def- 
inition and evolution of the perturbations, which is not our main concern. We have not 
attempted a full survey of this part of the literature, but have just cited useful papers that 
we happen to be aware of. By contrast, the cosmology literature on stochastic properties is 
very small for the case r^o < 1, and as we have mentioned it is out of touch with the relevant 
pure mathematics literature where the theory of random fields is discussed. The first serious 
treatment of stochastic properties is by Wilson in 1983 |14|. He developed the theory from 
scratch, and not surprisingly included only the sub-curvature modes which he knew were 
sufficient for the description of the non-stochastic properties. His notation is defective and 
much is left unsaid, but subsequent papers have not made basic advances in the formulation 
of the subject, though they have gone much further in calculating the cmb multipoles and 
comparing them with observation. We believe our referencing to be reasonable complete, 
as far as the cosmology literature on the stochastic properties is concerned. 

The layout of this paper is as follows. In Section 2 some basic formulas are given for the 
Robertson- Walker universe with 0, < 1. In Section 3 the standard procedure is described, 
and in the next section it is extended to the super-curvature modes. Inflation is discussed in 
Section 5, and the cmb anisotropy is treated in Section 6. In an Appendix we give various 
mathematical results in the sort of language that is familiar to us as physicists. 



2 Distance scales 



Ignoring perturbations, the universe is homogeneous and isotropic. There is a universal 
scale factor a{t), with t the universal time measured by the synchronized clocks of comoving 
observers, and the distance between any two such observers is proportional to a. 

According to the Einstein field equation, the time dependence of a is governed by the 
Friedmann equation which may be written 

Here i^T is a constant, H = a/a is the Hubble parameter, and O is the energy density 
measured in units of the critical density SH"^ /8itG. (As usual we set c = 1, and regard any 
cosmological constant as a contribution to the energy density as opposed to a modification 
of the Einstein field equation.) The spatial curvature scalar is 

= 6K/a^ (2) 

The distance defines the curvature scale; on much smaller scales space is practically 

flat, whereas on much bigger scales the effect of curvature is very important. From Eq. (|l|) 
the Hubble distance 1/H is a fraction (1 — f])^/'^ < 1 of the curvature scale. Even in the 
extreme case = 0.1 this makes the curvature scale about three times the present Hubble 
distance. 

We will set K = —1 so that a is the curvature scale. Then the case = 1 corresponds 
to the limit a oo, with physical distances like remaining constant. Note that the 
effect of curvature in a comoving region becomes neither more nor less important with the 
passage of time, since the curvature scale expands with the universe. 

We are concerned with the comoving region which is now the observable universe, 
bounded by the surface of last scattering of radiation emitted at very high redshift. This is 
close to the particle horizon of a matter dominated cosmology, unless there is a cosmological 
constant or some other non-standard contribution to the energy density. The coordinate 
distance of this particle horizon [|15| (ie., its distance in units of the curvature scale) is 
where sinh^ ^rph = U^^^ — 1. Even the smallest conceivable value JIq — 0.1 gives = 3.6, 
so effect of curvature is negligible except on scales comparable with the size of the observable 
universe. 

From Eq. the physical distance of the particle horizon is 



:i - J]o)"'/'^o"Vph (3) 



For r^o = 1 it is 2//q ^, and even for f^o = 0.1 it is only 3.8Hq ^. Thus it is not very much 
bigger than the Hubble distance H^^. 



3 Sub-curvature modes 



We are concerned with the first order treatment of cosmological perturbations. To this order, 
the perturbations 'live' in unperturbed spacetime, because the distortion of the spacetime 
geometry is itself a perturbation. 

The perturbations satisfy linear partial differential equations, in which derivatives with 
respect to comoving coordinates occur only through the Laplacian. When the perturbations 
are expanded in eigenfunctions of the Laplacian with eigenvalues —{k/a)'^, each mode (term 
in the expansion) decouples. 

Denoting the eigenvalue by — (fc/a)^, it is known ^ that the modes with real > 1 
provide a complete orthonormal basis for functions, and the usual procedure is to keep 



only them. Since they all vary appreciably on scales less than the curvature scale a we will 
call them sub-curvature modes. It will be useful to define the quantity 



g2 = A:^ - 1 (4) 

3.1 The spherical expansion 

Spherical coordinates are defined by the line element 

d/^ = a'^ldr'^ + sinh^ r{de^ + sin^ 9d(p'^)] (5) 

In the region r <^ 1 curvature is negligible and this becomes the flat-space line element 
written in spherical polar coordinates. The volume element between adjacent spheres is 
47rsinh^rdr, so for r ^ 1 the volume V and area ^ of a sphere are related by V = A/2. 
In contrast with the flat-space case this relation is independent of r, because most of the 
volume of a very large sphere is near its surface. 

Since the spherical harmonics are a complete set on the sphere, any eigenfunction 
can be expanded in terms of them. The radial functions depend only on r, and they 
satisfy a second order differential equation. As in the flat-space case, only one of the two 
solutions is well behaved at the origin, so the radial functions are completely determined up 
to normalisation. The mode expansion of a generic perturbation / is therefore of the form 

f{r,e,(j),t)=\ dqy^fkira{t)Zkim{r,0,4>) (6) 

•^0 Im 

where 

Zki^ = Iiki{r)Yi^{e,4>) (7) 
A compact expression for the radial functions is H, ||, g |l|] 



T{l + l+iq) i-l 



2 



which corresponds to the normalisation 

/ Hki {r)Uk', (r) sinh^ rdr = 6{q - q')5u' (9) 
JO 

The corresponding normalisation of the eigenfunctions is 

j ZklmZk'l'm'dV = 6{q - q')6ii'6mm' (10) 

where 

dV = sinh^ r sin OdrdOdip (11) 

is the volume element. 

As it stands Eq. (|8|) has a constant nonzero phase. It is convenient to drop this phase 
so that the function is real, and one then has the explicit expressions |[l^P| 

^ki = Nkitlki (12) 
tiki = g-^(sinhr)' (^^^\ cosiqr) (13) 



[2 , 

Nkl = \ -q^ 
V vr 



\sinhr dr , 

I 1 -1/2 



Ln=0 



(14) 



These expressions correct some misprints in [|l9, 



The un- normalised radial functions 11^; satisfy a recurrence relation |20| 



and the first three functions are 

f^k2 



1 

sinh r 
1 

sinhr 
1 



(/ + 1)2 + g2] + (2/ + 3) coth rUk,i+i 
sin(gr)" 



cos(gr) + cothr 



sm{qr) 



sinhr 



3 coth r cos{qr) + (3 coth^ r — — 1 



sm{qr) 



q 



(15) 

(16) 
(17) 
(18) 



The case Q = 1 corresponds to g — > oo with qr fixed, and in that limit Ilki(r) reduces 
to the familiar radial function, 



qjiiqr) ■ 



(19) 



Near the origin Ilki{r) has the same behaviour as ji{qr), namely 11^/ oc r', which ensures 
that the Laplacian is well defined there. The other linearly independent solution of the 
radial equation, which corresponds to the substitution cos{qr) —>■ sin(g'r) in Eq. (|^), has 
the same behaviour as the other Bessel function hi(qr) and is therefore excluded. 



3.2 Stochastic properties 

We are interested in the stochastic properties of the perturbations, at fixed time. To define 
them we will take the approach of considering an ensemble of universes of which ours is 
supposed to be one. 

The stochastic properties of a generic perturbation f{r,9,(j)) are defined by the set of 
probability distribution functions, relating to the outcome of a simultaneous measurement of 
a perturbation at a given set of points. From the probability distributions one can calculate 
ensemble expectation values, such as the correlation function for a pair of points ri,Oi,(j)i 
and r2,62,4>2, 

Cf^{fin,0i,cPi),f{r2,e2A2)) (20) 

and the mean square {f'^{r, 9, (/))). 

If the probability distributions depend only on the geodesic distances between the points, 
the perturbation is said to be homogeneous with respect to the group of transformations 
that preserve this distance. (For flat space this is the group of translations and rotations, 
and for homogeneous negatively curved space it is isomorphic to the Lorentz group pl| .) 
Then the correlation function depends only on the distance between the points, and the 
mean square is just a number. 

Cosmological perturbations are assumed to be homogeneous, and except for the curva- 
ture perturbation that we discuss in Section 6 their correlation functions are supposed to 
be very small beyond some maximum distance, called the correlation length. 



An ergodic universe? 

If there is a finite correlation length, one ought to be able to dispense with the concept 
of an ensemble of universes, in favour of the concept of sampling our own universe at 
different locations. In this approach one defines the probability distribution for simultaneous 
measurements at N points with by considering random locations of these points, subject 
to the condition that the distances between them are fixed. The correlation function is 



defined by averaging over all pairs of points a given distance apart, and the mean square is 
the spatial average of the square. For a Gaussian perturbation in flat space this 'ergodic' 



property can be proved under weak conditions |22] and there is no reason to think that 
spatial curvature causes any problem though we are not aware of any literature on the 
subject. 

For the ergodic viewpoint to be useful, the observable in question has to be measured 
in a region that is big compared with the correlation length. This is the case for the 
distributions and peculiar velocities of galaxies and clusters, where surveys have been done 
out to several hundred Mpc to be compared with a correlation length of order 10 Mpc, and 



accordingly the ergodic viewpoint is always adopted there |23]. However, even a distance 
of a few hundred Mpc is only ten percent or so of the Hubble distance Hq^, and therefore 
at most a few percent of the curvature scale (1 — flo)~^^'^HQ^ . Thus galaxy and cluster 
surveys do not probe spatial curvature. The only observables that do, which are the low 
multipoles of the cmb anisotropy, are measured only at our position so there is no practical 
advantage in going beyond the concept of the ensemble even if the mathematics turns out 
to be straightforward. 

In addition to the interpretation that the ensemble corresponds to different locations 
within the smooth patch of the universe that we inhabit, there are two other possibilities. 
One is that the ensemble corresponds to different smooth patches, which are indeed supposed 



to exist both in 'chaotic' |24] and bubble nucleation |25, 26, 27, 28 1 scenarios of inflation. 
The other, adopting the usual language of quantum mechanics, is to regard the ensemble 
as the set of all possible outcomes of a 'measurement' performed on a given state vector. 
A concrete realization of this 'quantum cosmology' viewpoint is provided by the hypothesis 
that the perturbations originate as a vacuum fluctuation of the inflaton field, which we 
consider later. 

3.3 Gaussian perturbations 

It is generally assumed that cosmological perturbations are Gaussian, in the regime where 
they are evolving linearly. A Gaussian perturbation is normally defined as one whose 
probability distribution functions are multivariate Gaussians [29, 22, 30 1, and its stochas- 



tic properties are completely determined by its correlation function. The perturbation is 
homogeneous if the correlation function depends only on the distance between the points. 

The simplest Gaussian perturbation is just a coefficient times a given function, the co- 
efficient having a Gaussian probability distribution. A more general Gaussian perturbation 
is a linear superposition of functions l2^ 



fir,e,(^)=Y,fnXnir,e,cP) (21) 

n 

with each coefficient having an independent Gaussian distribution. Its stochastic properties 
are completely determined by the mean squares (/^) of the coefficients. (For the moment 
we are taking the expansion functions Xn to be real, and to be labelled by a discrete index.) 
The correlation function corresponding to the above expansion is 

{f{n,9i,4>i)f{r2,e2,^2)) = Y.ifn)Xn{ri,ei,4>i)Xn{r2,02A2) (22) 

n 

For it to depend only on the distance between the points requires very special choices of 
the expansion functions, and of the mean squares 

It is very important to realise that the functions in such an expansion need not be linearly 
independent. Suppose for example that X^ = Xi + X2, and that (/I) is much bigger than 
(fi) and (/I)- Then most members of the ensemble are of the form / =constX3, which 
would clearly not have been the case if the function X3 had been dropped because of its 
linear dependence. 



So far all our considerations have been at a fixed time. The time dependence is trivial 
if we expand in eigenfunctions of the Laplacian, because each coefficient /„ then evolves 
independently of the others. Let us therefore replace the discrete, real expansion above 
by the complex, partially continuous expansion Eq. (^). The coefficients now satisfy the 
reality condition /^^^ = fki-rm and a Gaussian perturbation is constructed by assigning 
independent Gaussian probability distributions to the real and imaginary parts of the co- 
efficients with m > 0. We demonstrate in the Appendix that the correlation function being 
dependent only on the distance between the points is equivalent to the mean squares of 
their real and imaginary parts being equal, and independent of I and m. One can therefore 
define the spectrum of a generic perturbation / by 0] 



27r2 



(23) 



The q dependence of this expression has been chosen to give the simple form Eq. (30) for 
the mean square perturbation. In the flat-space limit (7 ^ oo it reduces to 

ifklmfk'l'm') = -^'Pf{k)S{k - k')5iii5^m' (24) 
The correlation function is given by 

( 2 ^i\ ^f^^)H^klm{n,Oi,(t>i)Zklm{r2,02A2) (25) 
^ ' Im 

Taking one of the points to be at the origin, only the term I = m = survives, and using 

Zkooir,e,(P) = {27r^)-^/^qsm{qr)/{qsmhr) (26) 

one finds 

Note the appearance in this expression of the logarithmic interval 

f = TT^ (28) 
k 1 + 

The flat-space limit is g ^ 00 with qr fixed, and the correlation function then reduces to 

, N 1"°° „ ,,^s'm(kr)dk 

(There is no distinction between k and q in the flat-space limit, and whenever we consider 
that limit we will use the the symbol k.) Setting r = gives the mean square value 

/CO At. 
-Vfik) (30) 



The flat-space limit is 



poo 

cm ^ if) = i ^Vf{k) (31) 







Since the spectrum is positive and | sin(qr)| < qr, the flat-space correlation function is 
never bigger than its value at r = 0, but without knowing the spectrum one cannot say 
more. The situation is very different for the curved-space expression Eq. (p7|), because it 
contains an extra factor r/sinhr. It follows that by expanding a perturbation in terms of 
sub-curvature modes one obtains a correlation function bounded by 

^ < ^— (32) 
^/(O) sinhr 



In order for (/^) to be well defined, the spectrum must have appropriate behaviour at 
q = oo and 0. As g — > oo one needs V ^ 0. As g ^ one needs 7^ ^ in the flat case, but 
only q^Vf{k) — > in the curved case. 

Note that in the curved case the limit g — > does not correspond to infinite large scales, 
but rather to scales of order the curvature scale. This means that one cannot tolerate a 
divergent behaviour there (unless of course the curvature scale happens to be larger than 
any relevant scale, in which case we are back to flat space). 

For future reference, we note that most other authors have used a different definition of 
the spectrum. This is usually denoted by Pj, and it is related to our Vf by 



With this definition. 



^/(^) = '^-^^^Pfik) (33) 



4 Including the super-curvature modes 

In the last section we found that the usual procedure, which includes only sub-curvature 
modes, generates a Gaussian perturbation whose correlation function necessarily falls off 
faster than r/ sinhr. This reflects the fact that each super-curvature mode varies strongly 
on a scale no bigger than the curvature scale. A random superposition of such modes will 



hardly ever be nearly constant on a scale much bigger than the curvature scale [31|, which 
is precisely what the lack of correlation on large scales is telling us. 

Faced with this situation one can might think that the lack of correlation on super- 
curvature scales is just a mathematical fact, inherent in the nature of homogeneous neg- 
atively curved space. Certainly it is not trivial to construct a function exhibiting super- 
curvature correlations. Consider, for instance, the following construction; throw down ran- 
domly spheres with radius much bigger than the curvature scale and fill them uniformly 
with galaxies, leaving the rest of the universe empty. Then one might think that a typical 
observer will see a uniform distribution of galaxies out to a distance much bigger than the 
curvature scale, making the correlation function almost flat out to such a distance. But this 
is incorrect, because according to the line element Eq. (^) most of the volume of a sphere 
is near its edge, and so is a typical observer.^ 

This example notwithstanding, correlation on arbitrarily large scales is possible, and is 
achieved simply by including the super-curvature modes. 

For — 1 < < the analytic continuation of the radial function 11^/ is purely imaginary, 
and for convenience we drop the i factor. Thus the super-curvature modes are deflned by 

liki = Nkitiki (35) 

/ — 1 d \ '"""^ 

nfc, = |gr2(sinhr)' (^— cosh(|g|r) (36) 

= \Mq\ (37) 

V TT 

,-1/2 

{I > 0) (38) 



,'2 , , 



Ln=l 

The recurrence relation Eq. ( [l5| ) is still satisfied, and the first three functions are 

'sinh{\q\r 



sinhr 



\Q\ 



(39) 



One of us (DHL) is indebted to R. Gott and P. J. E. Peebles for pointing out this fact. 



nfc2 



sinhr 
1 

sinhr 



, ,, , , , sinhflglr) 
cosh(|(;|rj + cotnr 



1 2 2 smh(|g|r) 
-3cothrcosh(|g|r) + (3coth r — q — 1) j— j 



(40) 
(41) 



At large r the super-curvature modes go like exp[— (1 — |g|)r]. Because the volume element 
is dV = sinh^ r sin ddrdOdip the integral over all space of a product of any two of them 
diverges. As a result they are not orthogonal in the sense of Eq. (|9|), let alone orthonormal. 
In any finite region of space (and of course we are only going to do physics in such a region) 
they are not even linearly independent of the sub-curvature eigenf unctions, since the latter 
are complete (for the set of functions defined over all space). None of this matters for 
the purpose of generating a Gaussian perturbation. 

The super-curvature modes add an additional term to the expansion Eq. 



e, (I)) = f' d{iq) J2 fklmZkimir, e, (I)) (42) 



Im 

Let us define the corresponding spectrum by analogy with Eq. ( p3| 

2^2 



ifklmfk' 



I'm') = |g|(g2 + l) ^/(fe)'^(lgl - WD^ll'^mm' ("1 < < 0) (43) 



We show in the Appendix that the correlation function remains well defined, and dependent 
only on the distance between the two points. Taking one of them to be at the origin only 
the I = mode survives and the super-curvature contribution to the correlation function is 
seen to be 

tSC. N _ d/c sinh(|g|r) 







The super-curvature contribution to the mean square is 







Unified expressions including all modes 

The use of q in the mode expansion Eq. (|6|) is natural for the sub-curvature modes, and we 
are using in this paper to facilitate comparison with existing literature. Unified expressions 
including all modes on an equal footing would use k in the mode expansion, so defining new 
coefficients fkim- One would then have the following expressions, which include both sub- 
and super-curvature modes. 

/(r,0,0,t) = / dkJ2fkUt)Zkim{r,e,^) (46) 



Im 

k\q 

= (48) 

Jo k gsmhr 



ifklmfk' I'm') = jr^,Vf{k)5{k-k')5u'6^m' (47) 



4.1 Very large super-curvature scales 

The contribution to the correlation function from a mode with k^ = 1 + q^ is 



if{r) oc exp(— fc^r) 



(49) 



for r ^ 1. Thus the correlation length, in units of the curvature scale a, is of order k~'^. This 
is in contrast with the flat-space case, where the contribution from a mode with /c <C 1 gives 
a correlation length of order 1/k. The difference can be understood in terms of the different 
behaviour of the volume element, in the following way. In both cases, the r dependence is 
that of the / = mode, and as long as r is small enough that the mode is approximately 
constant the divergence theorem gives 

r d/ ^ -fcVV(r) ^^^^ 
/ dr A{r) 

where V is the volume within a sphere of radius r and A is the area of this sphere. In flat 
space the right hand side is equal to — (A;r)^/3, so it is small out to a distance r ~ l//c. In 
curved space, when r ^ 1, the line element Eq. (^) shows that most of the volume of the 
sphere is near its edge, and the right hand side becomes equal to —k'^r/2, which is small 
out to a distance r ~ /c~^. 

In addition to being significant in its own right, the correlation function determines 
other physically significant quantities. One, which is relevant for quantities like the density 
perturbation, is the o"j(r), the mean square of / after it has been smeared over a sphere of 
radius r, which is given by 

a}{r)=V{r)-^ J dVi J dV2e/(ri, ^i, 0i, rs, ^2, </>2) (51) 

where the integrations are within the spheres ri < r and r2 < r (we have taken advantage 
of homogeneity to evaluate locate the sphere at the origin) . Because most of the volume of 
the sphere is near it's edge, the behaviour Eq. ( |49| ) of the correlation function leads to the 
same behaviour for crj(r), so it too remains constant out to a distance r ~ Another 
quantity, which is relevant for the curvature perturbation TZ that we shall discuss in Section 
6, is the mean square after smearing over a geodesic surface of radius r. 

Airy I dAi I dA2^Ti{ri, 01, cl)i,r2, 02,(1)2) (52) 

On super-curvature scales it is a measure of the fractional perturbation in the curvature of 
the surface. One sees that it too shares the behaviour Eq. (|4^), and so is constant out to a 
distance r ~ k~'^. 



5 The cmb anisotropy 



Even if Qq < 1) spatial curvature is negligible on scales that are small compared with the 
Hubble distance. As a result, the only observational data that can be sensitive to curvature, 
even if JIq < 1; are the lowest few multipoles of the cmb anisotropy. Papers discussing the 
effect of curvature on these multipoles O, ^, 33, 34] appeared sporadically before 
they were measured by the COBE satellite [^j, and many have appeared since [|3^, 37, 
38] 



4C, 41 1 . All of these papers keep only sub-curvature modes. Here we consider both 



sub- and super-curvature modes. 
The multipoles are defined by 



AT(e) 
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(53) 



The dipole term w.e is well measured, and is the Doppler shift caused by our velocity w 
relative to the rest frame of the cmb. Unless otherwise stated, AT will denote only the 
intrinsic, non-dipole contribution from now on. 



If the perturbations in the universe are Gaussian, the real and imaginary part of each 
multipole will have an independent Gaussian probability distribution (subject to the con- 
dition a^j^ = a^-m)- The expectation values of the squares of the real and imaginary parts 
are equal so one need only consider their sum, 



Ci = {M') . (54) 

Rotational invariance is equivalent to the independence of this expression on m. 

Even if it can be identified with an average over observer positions, the expectation 
value Ci cannot be measured. Given a theoretical prediction for Ci, the best guess for 
\aim\'^ measured at our position is that it is equal to C;, but one can also calculate the 
variance of this guess, which is called the cosmic variance. Since the real and imaginary 
part of each multipole has an independent Gaussian distribution the cosmic variance of 
J2m I^ZmP is only 2/(2/ + 1) times its expected value, and by taking the average over several 
Vs one can reduce the cosmic variance even further. Nevertheless, for the low multipoles 
that are sensitive to curvature it represents a serious limitation on our ability to distinguish 
between different hypotheses about the C;. Any hypothesis can be made consistent with 
observation by supposing that the region around us is sufficiently atypical. 

The surface of last scattering of the cmb is practically at the particle horizon, whose 
coordinate distance is r/o with sinh^ 770/2 = ^ — 1- An angle 9 subtends at this surface a 



coordinate distance d given by |23] 



6 = ^{l-nor^^^^od=^{aoHonod) (55) 
Spatial curvature is negligible when d <C 1, corresponding to 

e < 30(1 - no)-^^'^no degrees (56) 
A structure with angular size radians is dominated by multipoles with 

I ~ 1/9 (57) 
one expects that spatial curvature will be negligible for the multipoles 



(58) 

This is the regime / ^ 20 if JIq = 0.1, and the regime / ^ 6 if = 0.3. 

This restriction need not apply to super-curvature modes with /c^ <C 1 because the 
spatial gradient involved is then small in units of the curvature scale. The contribution of 
these modes is called the Grishchuk-Zeldovich effect, and we discuss it later. 

The linear scale probed by the multipoles decreases as / increases, and for / ~ 1000 it 
becomes of order 100 Mpc. On these scales one can observe the distribution and motion 
of galaxies and clusters in the region around us. On the supposition that they all have a 
common origin, the cmb anisotropy and the motion and distribution of galaxies and clusters 
are collectively termed 'large scale structure'. 

A promising model of large scale structure is that it originates as an adiabatic density 



perturbation, or equivalently |42, 43, 44, 45| as a perturbation in the curvature of the 



hypersurfaces orthogonal to the comoving worldlines. This model has has been widely 



investigated for the case = 1 [46|, and recently it has been advocated also for the case 



< 1 [35, 39, 37|. In this paper we consider the model only in relation to the cmb 



anisotropy since the galaxy and cluster data are insensitive to spatial curvature. We note 



though that the full data set may impose a significant lower bound on fig |47]. 



5.1 The curvature perturbation 



The curvature perturbation is conveniently characterised by a quantity TZ, which is defined 
in terms of the perturbation in the curvature scalar by^ 

4(fc2 + 3)7^fc,„/a2 = 5i^il (59) 

In the limit — > 1, 

Ak'TZuiJa' = 6R^l (60) 

On cosmologically interesting scales, TZkim is expected to be practically constant in the 
early universe. To be precise, it is practically constant on scales far outside the horizon in 
the regime where Q{t) is close to 1 (assuming that the density perturbation is adiabatic) 



[43, 44, 45, 52]. During matter domination the former condition can be dropped, so that 
T^kim is constant on all scales until 0, breaks away from 1. After that it has the time 
dependence TZkim = FlZkim where TZkim is the early time constant value and 

^ ^ sinh^ r/ — St/ sinh r/ + 4 cosh ?7 — 4 
(coshr/ — 1)2 

with 

T] = 2{aHy^ = 2(1 - n)^/^ (62) 

Unless otherwise indicated, TZ will indicate the primordial value TZ from now on. 

During matter domination and before O breaks away from 1, the density contrast is 
given by 

/ a'^H'^ \ 5pkim 2, 



For = 1 this reduces to 



fc2 + 3y p 5 
, P 5 



TZkim (63) 



T^kim (64) 



In these expressions the density perturbation is evaluated on comoving hypersurfaces (it is 
often referred to as the 'gauge invariant' density perturbation). In the matter dominated 
era where they hold this is the same as the density perturbation in the 'synchronous gauge' 



with the 'gauge mode' dropped [p3 |. 



5.2 The Sachs- Wolfe effect 

Horizon entry occurs long after matter domination on scales 

ao/k > 20(l^o/i^)~^ Mpc (65) 

where h is the value of Hq in units of lOOkmsec"^ Mpc~^. As a fraction of the Hubble 
distance these are the scales 

k Qoh 



It follows from Eq. (|55|) that they correspond to multipoles 



/ < 300/1 (67) 



^The quantity TZ was called (pm by Bardeen who first considered it , TZm by Kodama and Sasaki [ [48| . 
It is equal to 3/2 times the quantity 6K/k^ of Lyth which is in turn equal to the of Mukhanov, 

Feldman and Brandenberger [^. After matter domination it is equal to —(3/5)$, where "I" is the peculiar 
gravitational potential (and one of the 'gauge invariant' variables introduced in |k2|). On scales far outside 
the horizon, in the case f2 = 1, it is the of |Q, and three times the of jsi]] . 



Assuming an initial adiabatic perturbation, these multipoles are dominated by the effect of 
the distortion of the spacetime metric between us and the surface of last scattering, which 
is called Sachs- Wolfe effect. If = 1 the Sachs- Wolfe approximation accounts for about 
90% of C/ at / = 10, and about 50% at / = 30 §. 

The Sachs- Wolfe effect is determined by the curvature perturbation. In the case ^Iq = 1 



it is given by [23, 46, 



AT{e)/T = -\n{me) (68) 
5 

where tjq = 2(ao^^o)^^ is the coordinate distance of the edge of the observable universe, 
taken to correspond to the surface of last scattering. (The right hand side is usual given as 
^$(r/oe) where ^ is the peculiar gravitational potential.) Using Eq. (19) the multipoles are 
therefore given by 

1 f°° [2 
aim = -T dkTZkimX -kjiirjok) . (^9) 
5 Jo V vr 

Using Eq. (^4|), the mean square multipoles Ci are therefore given by 

4.77- roo 

Ci = ^l ^3f{mk)Vn{k) (70) 

For < 1, keeping for the moment only sub-curvature modes, the Sachs- Wolfe effect 
is given by [|, || |3| 

aim = - / dqUkimqlki (71) 
1 6 n° 

qlki = -Um{vo) + t drUMir)F'{ijo-r) (72) 
5 b Jo 

Here TZkim is evaluated well before O breaks away from 1 (when it is constant), F is given by 
Eq. ( |6l| ) and r]Q = 2(1 — rio)^''^ is again the coordinate distance of the edge of the observable 
universe. The mean square multipoles are given by 

/oo Ah 
-jVn{k)ll (73) 

When k ^ 1, I^i tends to a finite and nonzero limit for each /. This means that the C; are 
finite provided that q^V-ji 0. 

Two things should be noted about the regime g ^ in the curved space case. First, 
all multipoles receive contributions from this regime; in contrast with the flat case the 
quadrupole does not dominate as is claimed in pTI). Second, the limit q corresponds 



to scales of order the curvature, not to inflnitely large scales as is claimed in |^, 41]. Because 
of this last fact, one cannot tolerate a divergence of the C/ as g — > (unless the curvature 
scale is much bigger than any scale of interest in which case one is back to the fiat-space 
case) . 

5.3 Inflation and horizon exit 

It is widely supposed that the hot big bang is preceded by an era of infiation, during which 
gravity is by definition repulsive. A very attractive hypothesis is that the curvature per- 
turbation originates as a vacuum fluctuation during infiation, so that the ensemble average 
appearing in the definition of the spectrum (Eq. (p3|) ) is just the vacuum expectation value. 
Made originally for the case = 1 |55, IsTI, HI, [5^, 57 1, this hypothesis was later extended 



to the case r^o < 1 by Lyth and Stewart Before discussing it, let us see how infiation 
works with special reference to the case J^o < 1- 

The Hubble distance is usually termed the horizon (to be distinguished from the 
particle horizon), and the comoving length scale a/k associated with a given mode is said 



to be outside the horizon if aH/k > 1, and inside the horizon if aH/k < 1. The evolution 
of perturbations outside the horizon is very simple, because it is not affected by causal 
processes. Instead, the perturbation evolves independently in each comoving region p6| . 

Super-curvature scales, a/k > 1, are always outside the horizon (from Eq. but sub- 
curvature scales can be either outside or inside it. In the usual cosmology where gravity 
is attractive, aH = a decreases with time and at each epoch some scale is entering the 
horizon. The Hubble scale H^^ is entering the horizon now, and smaller scales entered the 
horizon earlier. Also, from Eq. (|^), Q is driven away from 1 as time passes, so that |1 — r2| 
must have been extraordinarily small at early times even if it is not small now. 

Inflation may be defined as an early era of repulsive gravity, when aH = a increases 
with time, and it is widely supposed that such an era preceded the hot big bang. At 
each epoch during inflation some scale is leaving the horizon, and as time goes by is 
driven towards 1. The standard assumption is that inflation occurs because the scalar field 
potential dominates the energy density, which falls slowly with time owing to the evolution 
of one of the scalar fields, termed the inflaton field. Constant energy density corresponds 
to (X and combining this dependence with Eq. (||) gives, for the case < 1, 

a = Hsinh{Ht) (74) 

and 

H = Hcoth{Ht) (75) 
After Q has been driven to 1, H achieves the almost constant value H, and 

a oc exp{Ht) (76) 

It is related to the scalar field potential V (in turn practically equal to the energy density) 

by 

= fmJ,fV (77) 

(After many Hubble times the Hubble constant could vary appreciably, in which case H 
denotes the value before this happens.) 

This evolution of the scale factor is modified if the energy density is rapidly decreasing 
with time; in particular, inflation might begin with a 'coasting epoch' during which 17 is 
almost constant Q |5^. We shall not consider that case. 

In the case f^o = !> inflation is usually held to solve at least three problems that arise if 
the hot big bang extends back to the Planck scale. Let us briefly recall them, and consider 
whether they are still solved if f^o < 1. For the moment we are discounting the bubble 
nucleation model of inflation. 

1. The harmful relic problem Without inflation it is difficult to avoid harmful relics of 
the early universe like monopoles or gravitinos. The possibility of avoiding them by 
inflation does not depend on the value of f^o- 

2. The flatness problem As we go back through time during the hot big bang, $7 is driven 
very close to 1 so that we have a fine tuning problem, the flatness problem. It is solved 
provided that inflation lasts long enough that 17 has been driven away from 1 again 
as we go back to the beginning of inflation. Prom Eq. (^) (or its 17 > 1 counterpart), 
the era at the beginning of inflation during which 17 is significantly different from 1 
has a typical duration of order where H is related to the scalar field potential by 
Eq. ([TTI). Prom Eq. (|l|), 17 has its present value 17o at the epoch during inflation when 
the observable universe (or to be precise, the comoving scale presently equal to the 
Hubble distance) leaves the horizon. Thus, 17o will be different from 1 if this epoch 
occurs near the beginning of inflation on the timescale H~^, but close to 1 otherwise. 



3. The homogeneity (horizon) problem Without inflation, the observable universe is far 
outside the horizon (Hubble distance) at early times. This means that causal processes 
cannot determine the initial conditions, which is usually held to be a problem, termed 
the 'horizon problem'. If is close to 1, the observable universe is typically far 
inside the horizon at the beginning of inflation, which solves the horizon problem and 
is usually said to 'explain' the homogeneity of the observable universe. Inflation with 
r^o < 1 cannot solve the horizon problem because the observable universe (or to be 
precise the comoving length presently equal to the Hubble distance) never occupies 
less than a fraction 1 — of the Hubble distance. 

However, no causal mechanism has ever been proposed for actually establishing homo- 
geneity at the beginning of inflation, even after the horizon problem has been solved. 
It seems to us therefore that the 'horizon problem' is a red herring, and that one 
should therefore look elsewhere for an explanation of the homogeneity of the universe. 
For the case r^o — 1 a fruitful avenue seems to be the following [^]. As smaller and 
smaller scales are considered one expects to find homogeneity below some minimum 
scale, but this is not the Hubble distance even though that is the only scale available 
at the classical level. Rather it is the scale, available only at the quantum level. 



57r/ \mpi, 

(we are setting h as well as c equal to 1). Indeed, within a volume with this radius, 
even the vacuum fluctuation of a massless scalar field generates energy density and 
pressure of order p^^'^, which would spoil inflation. As in the case of flat spacetime this 
vacuum contribution to the energy density is to be discounted (ie., one has to solve 
the cosmological constant problem by fiat at our present level of understanding). But 
one cannot allow a significant occupation number for the particle states defined on 
this vacuum. In other words, if f^o = 1 the universe has to be absolutely homogeneous 
at the classical level, on scales smaller than p"^/^. This guarantees the homogeneity 
of the observable universe at the classical level, provided that inflation starts at least 
\ln{mpi/Hi)] Hubble times before the observable universe leaves the horizon, where 
Hi is the value of H at this latter epoch. In order to respect the isotropy of the cmb 
one requires (Hi/mpi)^^'^ < |59|, and the bound is saturated in typical models 



of inflation. Thus, homogeneity of the observable universe is typically guaranteed if it 
leaves the horizon more than 7 or so Hubble times after the beginning of exponential 



inflation |46]. 

If f^o < 1 it is unclear how to define the vacuum as we discuss below, but with 
the mathematically simple conformal vacuum the vacuum fiuctuation again generates 
an energy density and pressure of order on the scale d. The criterion that this 
should not spoil the infiationary behaviour Eq. (^) is that d~^ be much less than 
the critical density, which requires as before d > {H/mpi)^/^H-^ . But now the 
observable universe is never far inside the Hubble distance H~^, so its homogeneity 
is not guaranteed by this type of argument. A different avenue would be to invoke 



quantum cosmology, along the lines of |60] which however deals only with the case 
> 1. 

The bubble nucleation model of inflation 

All of the above discussion assumes a classical evolution for the infiaton field, leading to a 
smooth evolution of S7(t). It might happen, however, that the scalar field potential allows 
quantum tunneling in scalar field space at some point during inflation. In that case a bubble 



of scalar fleld can form, whose interior is an J7 ^ 1 universe |25, 26, 28 1. Provided that 



the scalar field potential is still fiat enough, ft will again be driven to 1. 



If r^o turns out to be less than 1 the bubble nucleation model will be very attractive. 
Homogeneity is automatic. Also, CIq is determined by the form of the scalar field potential 



and can easily be less than 1 |2(:]. Assuming the usual 'chaotic' scenario for the beginning 
of inflation 62|, the inflaton field rolls slowly down a valley in scalar field space, and 

then the bubble nucleation model might correspond to sideways tunneling out of this valley 
The only problem would be to find a potential of the required form that looks sensible 



in the context of modern particle theory; as has recently been pointed out |58, this 
constraint makes it difficult even to find a potential that leads to ordinary, non-tunneling 
inflation. 

5.4 Sub-curvature contributions and the vacuum fluctuation 

During inflation, the curvature perturbation is related to the perturbation 6(j) of the inflaton 
field by [|7|, |, 1^ 

7^ = -{H/^)6(j) (79) 

where the dot denotes differentiation with respect to time t. This expression holds at all 
epochs, not just when TZ is constant. In it, 6<j) is defined [^] on hypersurfaces which have 
zero perturbation in their curvature scalar (it is often called the 'gauge invariant' inflaton 
field perturbation). 

A very attractive hypothesis is that 6(j) originates as a vacuum fluctuation, so that the 
ensemble average appearing in the deflnition of the spectrum (Eq. (^)) is just the vacuum 
expectation value |55, 43, 57, 35 1. 



The vacuum fluctuation during inflation also generates a spectrum of gravitational 
waves, which is well understood for the case = 1 |66|, and is under investigation for 
the case r^o < 1 P^. We will not consider it here. 



To calculate the vacuum fluctuation one uses quantum field theory in negatively curved 
space and the first step in setting up this theory is to expand in terms of the 

sub-curvature mode functions. In this context there is no question of including additional 
modes, because many results of quantum field theory (such as the vanishing of field com- 
mutators outside the light cone) depend essentially on the fact that one is using a complete 
orthonormal set. As a result the spectrum predicted by the vacuum fiuctuation will include 
only sub-curvature modes. 

The same restriction holds for the fluctuation in any quantum state that is homogeneous 
(with respect to the group of coordinate transformations leaving the distance between each 
pair of points invariant). But one can give the inflaton field perturbation any desired 
stochastic properties by choosing a suitable quantum state (pure or mixed), and in particular 
one can generate an arbitrary homogeneous Gaussian perturbation. The absence of super- 
curvature modes in the vacuum fluctuation prediction is not a feature of quantum fleld 
theory per se. 

The coefficients in the mode expansion of the quantum field (pkim satisfy the classical 
field equation (in the Heisenberg representation), which fixes them up to a one-parameter 
ambiguity once a convention is made for their normalisation. Breaking this ambiguity is 
equivalent to defining the vacuum. In the case = 1 each mode starts out well inside the 
horizon, where the spacetime curvature is negligible. In that case the vacuum is defined to 
be the usual flat spacetime vacuum, and assuming the usual slow roll conditions one flnds 
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In this expression V{(p) is the inflaton potential, and the right hand side is to be evaluated 
at the epoch of horizon exit k = aH. It gives an almost scale-independent result for typical 
models of inflation. 



In the case < 1, without bubble nucleation, it is not clear how to define the vacuum 
because a given scale is never far inside the horizon. The mathematically simplest choice 
is the 'conformal vacuum', and using it one finds (after suitably generalising the slow roll 
conditions) that Vtz is still given by the above expression [Q, (For the special case of 
a linear potential this result has been reproduced recently, using a different calculational 
technique p9[.pl) 

In the bubble nucleation model the quantum state of the inflaton field perturbation 
inside the bubble can be calculated [EtI, |26|, 1671, and it is found not to be in the conformal 



vacuum. As a result [g6|, VTi{k) is multiplied by a factor coth(7r(jf) compared with 
Eq. d^. 

Comparison of the vacuum fluctuation with observation 

If r^o = li Eq- (|70|) with a scale independent Vn gives 



1(1 + m = ^Vn (81) 

The prediction is that the left hand side is scale independent, which is consistent with 
observation |^. There is however room for considerable scale dependence; defining n by 
Viz q""^) the allowed range is 0.6 < 1.4. The magnitude is small, 

l{l + l)Ci = 8.05 X 10-10 Qrms-ps ^ g Q ^ ^Q-io (32) 

20 /xK 

which corresponds to 

~ 3 X IQ-^ (83) 

For r^o < 1 Eq. ( |73|) has to be evaluated numerically. With the flat spectrum coming 
from the conformal vacuum assumption, 1(1 + 1)C; has in general a negative slope for 
0.1 < Oq < 1 (for / in the range I < 10 where the Sachs- Wolfe effect dominates, and in 



which curvature can be significant) Q With present data the slope is not strong 



enough to rule out any value of though better data will probably rule out low values. 

As already noted, the spectrum is not flat in the bubble nucleation model, but rather 
is proportional to coth(7rg). It turns out however [28| that if JIq is substantially below 1 



the integral in Eq. ([72D dominates, with peaking at > 1 even for the quadrupole (and 
at higher values for higher multipoles). As a result the bubble nucleation prediction is not 
significantly different from the flat spectrum prediction, when cosmic variance is taken into 
account. 

Power law parameterization of the density perturbation spectrum 

The spectrum V-ji of the curvature perturbation is directly related to the Sachs- Wolfe effect, 
and from this viewpoint the assumption that it is effectively flat seems natural. This 
assumption was not, however, the assumption made in the literature before the (very recent) 
advent of the vacuum fluctuation prediction. Rather, it was assumed that the spectrum 



Ps (defined by Eq. (33)) of the density perturbation is proportional to q. This choice is 
equivalent to the flatness of V-ji for 0,q = 1, but otherwise it is equivalent to 

g^(l + g ^) 
(4 + g2 



-Pn oc '-r^^^ (84) 



^The authors of l39[ d o not establish the identity of their result with the earlier one, but it follows by 
evaluating Eq. (2) of ]37| (multiplied by 16n/m%i to bring the conventions of jsij into line with the usual 
ones) during matter domination before fl breaks away from unity. To do this one has to replace 1 — S7o by 
1 — n ^ 1 in the quantity Wi/ci, leading to Wi/ci ~+ {2/5)(aH)~^ . Remembering that the energy 
density scales like during matter domination and like a"** during radiation domination, one then indeed 
reproduces the spectrum of the energy density given by Eqs. ( |6^ ) and ([80[). 



The right hand side tends to 1 in the hmit g — > oo of neghgible curvature, but is much less 



than 1 for < 4. With this parameterization, l{l + l)Ci acquires a positive slope [36| for 
0.1 < J7o < 1, comparable in magnitude with the negative slope of the vacuum fluctuation 
prediction. 

The main cause of this difference is the factor + coming from the relation between 
the density and curvature perturbations, and a similar result would probably be obtained 
if Vs were used instead of Ps, or k instead of q. In other words, the predicted C/ can be 
regarded as coming simply from a linearly rising density perturbation spectrum, as opposed 
to a flat curvature perturbation spectrum. These two parameterizations are equivalent if 
ilo is equal to 1, but if f^o is significantly smaller than 1 the first one gives less power on 
small scales, leading to a significantly different prediction for the Q's. The central values 
of the present data points lie between the two predictions, and the present error bars are 
big enough that they are indistinguishable. But in the future the data should be able to 
distinguish between the two parameterizations, ruling out one or both of them for small 
values of ^Iq. 

5.5 Super-curvature scales and the Grishchuk-Zeldovich effect 

Like any statement in physics, the statement that the inflaton field is in the vacuum will be at 
best approximate, and its validity will presumably depend on the scale under consideration. 
When considering departures from it there is no reason to exclude super-curvature scales. 

The contribution of super-curvature scales to the mean square multipole C/ is just the 
extension of Eq. ( [73[ ) to super-curvature scales, 

CfC = fvn{k)ll (85) 

Requiring that the super-curvature contribution be no bigger than the total gives (Eq. (|8 

fCf" < 10-10 (86) 

The quantities /|; have not been calculated in the super-curvature regime < /c^ < 1, 
but as discussed below they are proportional to k'^ near k = 0, and we noted earlier that in 
the regime k^ > 1 they all peak at a value k'^ > 2. It therefore seems reasonably to suppose 
that for a fairly flat spectrum VTz{k), the super-curvature regime does not contribute much 



to the mean square multipoles Q. In that case Eq. (SC) will provide no significant constraint 
on a flat spectrum. In other words, it will be difficult to detect the cutoff below k^ = 1 that 
the vacuum fluctuation predicts. (The even more difficult task of finding an observational 



signal for the this cutoff without assuming that the spectrum is flat is discussed in [70|.) 

Now suppose, in contrast, that the spectrum rises sharply on some very large scale. Then 
there might be a big curvature perturbation, which would however have a big correlation 
length and so have a very small spatial gradient. If the correlation length is big enough, the 
gradient will be small enough to ensure that the curvature perturbation has no significant 
effect on the cmb anisotropy, even if it is quite is large. How big does the correlation length 
have to be for this to happen? 

For r^o = 1 this question was asked and essentially answered by Grishchuk and Zeldovich 
|10|. We will now briefly recall their argument, using the precise concept of the spectrum 



of the curvature perturbation in place of their more qualitative discussion. Then we will 
generalize it to the case < 1, and flnally discuss its physical significance in both cases. 

5.6 The Grishchuk-Zeldovich effect (Qq = 1) 

From Eq. (]3l|), a flat spectrum Vu gives a logarithmically divergent result for (7^^), but 
since V-jz ~ 10"^ one has to go to a huge scale to see any effect. Taking the comoving 



small-scale cutoff to be the Hubble distance at the end of inflation (which is equivalent to 
the usual procedure of dropping the contribution of the vacuum fluctuation to the energy 
density in flat spacetime), this scale as a multiple of the Hubble distance is 

~ exp(109 - 60) ~ exp(lO^) (87) 

If as expected the spectrum increases somewhat with scale this estimate will be sharply 
reduced, but it will still be a big number in typical models of inflation. 

A possible divergence of the geometry distortion associated with a nearly flat spectrum 



is interesting, because it suggests that the universe might be fractal on very large scales |61|. 
It cannot however be explored observationally, because if the spectrum is fairly flat, large 
scales give a negligible contribution to the cmb anisotropy. Here we are concerned with the 
quite different possibility, that the spectrum might rise sharply on some not-so-huge but 
still very large scale gvL- We therefore suppose that the spectrum has the form 



Such a contribution might originate from the vacuum fluctuation if the inflaton potential 
has a suitable form, but more plausibly it would arise because the vacuum assumption failed, 
or in other words because there were inflaton particles with momentum k < /cvL • 

Before calculating its effect, let us spell out the physical significance of this contribution. 
We will make the natural assumption that the homogeneous Gaussian perturbation under 
consideration exists in a patch around us, whose size is much bigger than the correlation 
length 

dvL = ao/kvL (89) 

The curvature perturbation is more or less constant in a region of size dvL- According to 
Gauss-Bonnet theorem, the distortion in geometry in this region (measured for instance by 
the departure of the sum of the angles in a geodesic triangle spanning it from 27r) is of order 
its cross-sectional area times the perturbation in the curvature scalar. From Eq. (|3l| ) this 
is of order 

4l^^^'^ ~ (^')'/' (90) 
In some patches the curvature perturbations will be positive, and in others negative. It 
does not make sense to consider a value (7^^) bigger than 1, because then the regions of 
positive curvature would close on themselves. 

The geometry distortion in the observable universe is of order UqSR^^^ , and so smaller by 
a factor feyL- As the correlation length increases with fixed (7^^), the geometry distortion in 
the observable universe decreases and so does the spatial gradient For both these reasons, 
the anisotropy caused by a given value of V-ji decreases. Let us calculate it. 

Since j/(x) ~ for small x, one sees from Eq. (|70|) that the quadrupole dominates. 
Using j2(x) = one finds 

The quadrupole measured by COBE is not significantly in excess of the typical values 
/(/ + l)Ci ~ 8 X 10"^*^ of the other multipoles (in fact it is somewhat smaller), so we 
conclude that is absent at this level. This means that 

dvL > 70(7^2)l/4//-l (92) 

As Grishchuk and Zeldovich pointed out, this bound on {TZ'^) becomes weaker as the 
scale increases. It can be of order 1 provided that (cf. [TljJ) 

^ > 70 (93) 



A bound similar to this one was already implied by upper limits on the quadrupole 
that existed two decades ago. But before COBE measured the actual values of the low 
multipoles there was always the possibility that the Grishchuk-Zeldovich effect might be 
present (ie., that the quadrupole might stick out above the other multipoles), indicating a 
big curvature perturbation on some very large scale. The somewhat disappointing fact that 
the Grishchuk-Zeldovich is absent seems not to have been noted anywhere in the copious 
literature on the COBE observations. 



5.7 The super-curvature scale Grishchuk-Zeldovich effect 

To generalize the Grishchuk-Zeldovich effect to JIq < 1 one needs to take spatial curvature 
into account, and to note that the limit of large scales corresponds to A; — > 0, not q ^ 0. 
This has not been done to date. The only relevant publications of which we are aware are 
1 71 1 where spatial curvature is ignored, and [p^, 41] where the appropriate limit is incorrectly 



assumed to be g ^ 0. In contrast with the case Qq = 1, the Grishchuk-Zeldovich effect for 
the case r^o < 1 could come only from a non-vacuum contribution. 

Consider therefore a contribution on very large super-curvature scales, and represent it 

by 

^6{lnk-ln /cvl) (94) 

We need to understand the physical significance of the corresponding curvature perturba- 
tion, which is different from the case 0,q = 1. From Section 4.1, the correlation length is 
now 

c?VL = ao/^VL (95) 

The curvature perturbation is practically constant in a region of this size. The perturbation 
in the geometry distortion of such a region is now 

4l5^('^ - k^tinY' (96) 

On the other hand there is now a distortion even in the absence of a perturbation, given by 
Eq. i), 

Thus (7^^)^/^ now measures the fractional perturbation in the geometry distortion, not 
the distortion itself. However the requirement that regions of space should not close on 
themselves is still equivalent to (Ji?)^/'^ < 1. 
An equivalent way of viewing 

(7^2)l/2 is that it measures the geometry distortion of 
a region with size oq. Since VLq is not extremely small this is the same as saying that it 
measures the geometry distortion of the observable universe. As the correlation length dvL 
is increased with {Ti?) constant, the geometry distortion of the observable universe does 
not decrease as it does in the = 1 case. One still expects, though, that for a given value 
of Vtz the effect on the cmb anisotropy will become smaller, because the spatial gradients 
become smaller. Let us see how to calculate it. 

As /c ^ 0, Ilfco 1, but the other radial functions are proportional to k. The normali- 
sation factor becomes 

kNki ^ Ni = \l^ (98) 
V vr 

kNki - Ni^J^[Y{{n''-l)\-^l^ {l>2) (99) 



n=2 



and 



^2 Y 

flfci (r) ^ — — ^ [sinh(2r) - 2r] (100) 
4 smh r 



The other radial functions follow from the recurrence relation Eq. (|l|). It is convenient to 
define 

Ui = lim Uki/k (101) 
Using these results, the contribution to the mean square multipoles becomes 

where 



1 ~ 6 f^'o ~ 

Bi = -ni{7]o) + - / Ui{r)F'{rjo - r)dr (103) 
5 5 Jo 



1 ^ . . 6 no ^ 
5 70 

When one increases the value of / under consideration, the scale above which these limits 
hold presumably becomes successively larger, so to actually calculate C^^ for a given value 
of fcvL one ought to use the full expression Eq. (85), but hopefully Eq. ( |102| ) will provide a 



reasonable estimate for small /. Since Ni and Bi are roughly of order 1 for low multipoles, 
it says very roughly that 

^ ^2^(7^2^ (104) 
The absence of at the level 10"^*^ therefore implies very roughly 

kyl> 10'' in') (105) 

Since (aoH)'^ = 1 — is supposed not to be tiny, this result is roughly 

dvL > 10^°(7^2)Fo-l (106) 

In words, the conclusion is that if the fractional geometry distortion is of order 1, its 
correlation length must be more than 10^^ Hubble distances. This result is not directly 
comparable with the Qq = 1 result, because it concerns the fractional, not the absolute, 
geometry distortion on the scale dvL- The quantity that measures the absolute geometry 
distortion is TZkim = kyiTZkim, and in terms of this quantity 

dvL >100(7^2) 1/5/7-1 (107) 

We see that if the absolute geometry distortion, in a region whose diameter is equal to the 
correlation length, is of order 1, then the correlation length must be at least two orders 
of magnitude bigger than the Hubble distance. This is essentially the same as the ^Iq = 1 
result. 

These estimates have been derived from the fact that the Grishchuk-Zeldovich effect 
cannot be bigger than the observed values of the multipoles. As we have not actually 
calculated the / dependence of we cannot say that the effect is definitely absent as 
in the = 1 case, because it might turn out that the dependence of Cvl mimics the 
dependence that of the data (roughly Ci oc It would be desirable to calculate the 

shape of C^^, both to check that this does not happen and to check the assumption that 
Bi is of order 1. 

In contrast with the case Qq = 1, the Grishchuk-Zeldovich effect is present in all of 
the low multipoles if J7o < 1- Indeed, it could even occur in multipoles / > 10 in which 
case the Sachs- Wolfe approximation would become inadequate to investigate it and a full 
calculation would be necessary. The necessary formalism to perform such a calculation is 
already in place |34], and it has already been used for the sub-curvature modes |3(:, 37, 69 1. 



The extension to the super-curvature modes raises no new issue of principle. 



5.8 The physical significance of the Grishchuk-Zeldovich effect 



In some of the literature |7l| , [4l[ , the absence of the Grishchuk-Zeldovich effect at a given 
level has been regarded as evidence that the smooth patch of the universe which we occupy 
extends beyond the edge of the observable universe. As we now explain, this is not the case. 



To make the simplest point first, it is clear from Eqs. ( [7l|) and (72) that the multipoles 
aim of the cmb anisotropy depends only on the curvature perturbation within the observable 
universe. This remains true when we consider their ensemble mean squares Ci. Strictly 
interpreted, the Grishchuk-Zeldovich effect just explores the effect of very small spatial 
gradients of the curvature perturbation, within the observable universe, on the hypothesis 
that the curvature perturbation is a typical realization of a homogeneous Gaussian random 
field. Recall that in this context 'homogeneous' means that the correlation function of 
the curvature perturbation depends only on the distance between two points, not on their 
location; but we are still talking about locations within the observable universe. 

If this hypothesis is indeed correct within the observable universe, one expects it to 
remain correct in some larger region. If this region is sufficiently big, one can introduce 
the concept of a correlation length as we did in the above discussion. By definition, the 
correlation function is more or less constant out to a distance of order the correlation length, 
only then falling off. Clearly 'sufficiently big' means bigger than the correlation length. But 
we will never know whether this picture is correct, because we will never know what lies 
beyond the edge of the observable universe (except by waiting for it to gradually recede, in 
comoving distance units). 

Homogeneity of the perturbation corresponds to the spectrum (defined by Eq. (^)) 
being independent of I and m. One can reasonably expect this property to fail when k 
becomes so small that the corresponding distance ao/ky^ becomes bigger than the size of 
the smooth patch of the universe around us, within which the perturbation is homogeneous. 
In that case the presence of the Grishchuk-Zeldovich effect on a given scale would suggest, 
though not really prove, that the smooth patch extends to the corresponding distance. But 
its absence says nothing. 

Beyond the smooth patch might be regions of the universe where the 'perturbations' 
become so big that it makes no sense to talk about a Robertson- Walker universe. If so 
the patch discussed in the last paragraph will have a periphery, within which the typical 
magnitude of the perturbations becomes bigger as one moves outwards (in contrast with the 
region within the patch, where the typical magnitude is by definition the same everywhere). 
The Grishchuk-Zeldovich effect tells us absolutely nothing about this periphery. 

For the case = 1 one might argue that information about the periphery is available, 
by taking the density perturbation to be the primary quantity rather than the curvature 
perturbation. From this viewpoint the large density perturbation in the periphery will gen- 
erate a large curvature perturbation in the observable universe, through the usual 'Coulomb 
law' solution of the Poisson equation Eq. (^), unless there is an accidental cancellation. 
However there does not seem to be any justification for it, and it does not work for < 1 
because according to Eq. ( |63| ) the Poisson equation does not hold. Rather, the density 
perturbation and the curvature perturbation become essentially the same on scales much 
bigger than the curvature scale. 

Inflation and the Grishchuk-Zeldovich effect 

A separate issue is whether the absence of the Grishchuk-Zeldovich effect tells us anything 
about inflation. This is clearly the case only if the scale /cvl can be related to inflation. 
Reference |71], which deals with the case = 1 (or at any rate ignores spatial curvature). 



accepts the usual dogma that 'the universe is smooth on some scale of order the Hubble 
distance at the beginning inflation'. Interpreting this to mean that /cvl ~ clH at the 



beginning of exponential (0 = 1) inflation, the absence of the Grishchuk-Zeldovich effect 
indeed tells us that inflation starts several Hubble times before the observable universe 
leaves the horizon. However, as discussed in Section 5.3 the usual dogma does not have any 
clear justification. 

6 Conclusion 

In this article we have drawn the attention of physicists to the incompleteness of the standard 
mode expansion for cosmological perturbations in an Qq < 1 universe. In order to generate 
the most general homogeneous Gaussian random field one should use Eq. (|47|) , which runs 
over all negative eigenvalues — A;^, whereas the standard expansion keeps only the modes 
with > 1. We have called these sub-curvature modes, because they vary appreciably 
over a distance less than the curvature scale, and we have called the modes with < A;^ < 1 
super-curvature modes. 

The fact that super-curvature modes are needed to generate the most general pertur- 
bation has been known to mathematicians for about half a century, so that their omission 
by cosmologists constitutes a remarkable failure of communication between the worlds of 
mathematics and science. This omission leads to perturbations which are practically un- 
correlated on scales bigger than the curvature scale. In contrast, a mode with A;^ <C 1 
corresponds to a perturbation with correlation length in units of the curvature scale. 

What nature has chosen to do is of course another question. For the case = 1 
the standard assumption is that the perturbations originate as the vacuum fluctuation of 
the inflaton field, and in 1990 this assumption was extended to the case ilo < 1 Q. The 
mode expansion of a quantum field runs only over sub-curvature modes, since they form a 
complete orthonormal set for square integrable function. As a result the vacuum fluctuation 
generates a perturbation which includes only these modes. 

Of course this is not a feature of quantum field theory per se, but of the assumption 
that the inflaton field is in the vacuum. This assumption might break down below some 
small value of k, corresponding to a correlation length much bigger than the curvature 
scale. We therefore ask whether a big perturbation with a very large correlation length 
could be detected through the cmb anisotropy. For the case ^Iq = 1, this question was 
asked and answered in 1978 by Grishchuk and Zeldovich ||l^, and here we have extended 
their discussion to the case J7o < 1 by including the super-curvature modes. We have 
given a formula for the cmb anisotropy due to a mode with k <^ 1, and have estimated 
its magnitude. By requiring that it be no bigger than the observed anisotropy we have 
estimated a lower limit on the correlation length for a perturbation of given magnitude. As 
in the case r^o = 1, the correlation length must be more than about two orders of magnitude 
bigger than the size of the observable universe, if the geometry distortion is of order 1 in a 
region whose size is equal to the correlation length. 

In contrast with the case f^o = 1, the Grishchuk-Zeldovich effect is present in all mul- 
tipoles, not just in the quadrupole. It would be interesting to evaluate its I dependence, if 
only to check that it does not mimic the observed dependence which is usually interpreted 
as coming entirely from the vacuum fluctuation. 

Appendix 

This appendix gives some mathematical results, in the sort of language that we as physicists 
are accustomed to. It deals with both the spherical expansion used in the text, and with 
an expansion using coordinates that slice space into flat surfaces which is more like the flat- 
space Fourier series. At the risk of being pedantic we give a rather full treatment, because 



even in the flat-space case there does not seem to be a reference that explains the basic 
concepts in a way that is accessible to most physicists. 

We refer the reader seeking a rigorous but more abstract treatment to [11|. 



The Fourier expansion 

In flat space the simplest approach is to use the Fourier expansion. In comoving coordinates 
r = (x, y, z) it is 

/(r) = (27r)-3/2 J d3q/q exp(iq.r) (108) 
The orthonormality relation is 

(27r)"3 J dV exp(-iq.r) exp(iq'.r) = ^^(q' - q) (109) 

where dV = d^r is the volume element. 

A Gaussian perturbation is obtained by assigning independent Gaussian probability 
distributions to the real and imaginary parts of the coefficients /q, but we need to en- 
sure translation and rotation invariance of the correlation function. Translation invariance 
is equivalent to the real and imaginary parts having the same distributions (same mean 
squares), for the following reason. Because exp[i(q.x)] exp[i(q'.x')] is a function of x — x' 
only if q = — q', translation invariance is equivalent to introducing a correlation between /q 
and /_q only, which because of the reality condition /* = /_q means a correlation between 
/q and /* only. This means that the phases of /q must be uncorrelated, which indeed means 
that the real and imaginary parts of /q must have the same distribution .[[] Let us therefore 
define the spectrum by 

27r2 

The correlation functions is then 

e/(r) = (27r)~327r2 / ^ exp(-iq.r)P/(g) (111) 
J (I 

Rotational invariance is clearly equivalent to the spectrum depending only on the magnitude 
of q, not its direction. Performing the angular integration one obtains Eq. (pS|). 

Using the well known expansion of a plane wave into spherical waves one can prove that 
the above definition of the spectrum is equivalent to the definition Eq. ( p^ in terms of the 
spherical expansion. The equivalence, and in particular the fact that the two definitions are 
the same except for the different delta functions, does not depend on the detailed form of 
the transformation between the spherical expansion and the Fourier expansion, but rather 
on the fact that it is unitary. 

The spherical expansion 

We first justify the claim made in the text, that the correlation function depends only on 
the distance between the points if the spectrum defined by Eq. ( p3| ) is independent of I and 
m. Let r, 9, (f) and r', 9', cj)' be the coordinates of a given point with respect to two different 
spherical coordinate systems. We saw earlier that the most general eigenfunction with 
eigenvalue —{k/af' is a linear combination of the functions Zkim{r,9,(f)). This is of course 
true in any coordinate system. Since Zkim{r, 0, 4>) and Zkim{r', 9', cj)') are both eigenfunctions 

direct way of seeing this is to work with the real form of the Fourier integral and note the identity 
cos a cos b + sin a sin b — cos(a — b) 



it follows that either of them can be expanded in terms of the other. Thus, there is a linear 
combination of the form 



Zklm{r,6,<j)) = ^ UijnVm'ZkVm'{r\0\(t)) (112) 
I'm' 



If the spectrum defined by Eqs. (|23D and (^) is independent of I and m, the correlation 
function given by Eq. (^) or its super-curvature analogue will be invariant under transfor- 
mations of the above form provided that the transformation matrix satisfies the unitarity 
property 

Ulml"m"{Ul'm'l"m")* = ^ll'^mm' (H^) 

r'm" 

Since the transformation takes a pair of points into arbitrary positions subject to the con- 
straint that the distance between them is fixed, the correlation function will then depend 
only on this distance. 

For the sub-curvature modes, unitarity follows from the fact that the transformation 
takes one orthonormal basis (for the subspace of eigenfunctions with a given eigenvalue) 
into another. Let us see this explicitly. Orthonormality for the whole space gives the 
the coefficients of the expansion as 

Siq - qWLi'm' = I Zh'm'ir', 9\ <P')Zkim{r, 6, 0)dV (114) 



where the primed coordinates are regarded as functions of the unprimed ones and the volume 
element is defined by Eq. (11). Now consider the inverse transformation. 



Zki'm'{r',0',cl)') =^Vy^,i^Zkim{r,0,(t)) (115) 

Im 

The coefficients are given by 

- <l)yum'lm = I Zhmir, 6, 0)Z,;,„,(^', <P'W (116) 



where now the unprimed coordinates are regarded as functions of the primed ones. But 
as the integration goes over all space one can just as well integrate over the unprimed 
coordinates and regarded the primed coordinates as the dependent ones. By comparing 
Eqs. (114) and ( |116| ) it follows that the transformation is indeed unitary. 



iyiml'm') — ^Vm'lm (H'^) 

This proof of unitarity does not work for the super-curvature regime, because we invoked 
orthonormality to obtain Eq. ( |114| ) for the matrix element Ui^y^i- There is however an 
alternative expression that remains well behaved in the super-curvature regime, obtained 
by substituting into Eq. ( |112| ) the definition of the Z's, and remembering that the spher- 
ical harmonics Yim ^ complete orthonormal set on the sphere. Choosing any sphere 
r' =constant this expression is 

Tikl'{r')Utmi'm' = J Uki{r)YiUO,m'm'{0',(l>') Sin e'de' del)' (118) 

where the unprimed coordinates are regarded as functions of the primed ones. (For sub- 
horizon modes the original expression Eq. ( |114| ) is recovered if we multiply both sides of 
Eq. (118) by Ilk'i'{r') sinh^ r'dr' and integrate over < r' < cxd, but continuing to imaginary 
q and q' causes the integrals on both sides to diverge at r' = oo.) Similarly, choosing any 
sphere r =constant one has 

n^K^V^/'m'^m = / Uki'ir')Y,m'{9', <P')Yiim sinededcj) (119) 



In contrast with the original expressions Eqs. (114) and ( |116| ), the radial coordinate on the 
right hand side of these new expressions runs only over a finite range. When we analytically 
continue to imaginary q, the radial functions pick up a factor i which cancels, so the equality 
iylrni'm'^ ~ ^Vm'im ^'^^^ established for real q remains valid. Note that this does not 
work for complex eigenvalues, indicating that they are not allowed. 

The above discussion is a generalisation of the familar demonstration that a rotation 
around the origin acts on 1/^ with a finite dimensional unitary matrix acting on the m index 
alone. Because of the finite dimensionality, this provides a rigorous proof that invariance 
under such rotations is equivalent to the independence of the spectrum on m. The extension 
to arbitrary rotations and translations involves an infinite sum over /, and as we have not 
discussed its convergence the discussion is not rigorous. Sasaki and Tanaka [p2| have recently 



demonstrated that the infinite sum over I that occurs in Eq. (25) is uniformly convergent, 
which makes the above derivation rigorous. 

The flat-surface expansion 



The flat-surface expansion uses coordinates defined by the line element 17, 14] 

dl^ = {a/zf{dx^ + dy^ + dz^) (120) 

The surfaces of constant z are flat. Any point in space can be chosen as the point x = y = 0, 
z = 1 and curvature is negligible in the region around that point |x| <C 1, \y\ <C 1, \z—l\ <C 1. 

Since a sphere of infinite radius is flat, one can think of the flat surfaces z =constant 
as spherical wave fronts, originating from a point at z = 0. Note that these surfaces have 
an 'inside' and an 'outside' even in the limit where they are flat, because geodesic surfaces 
such as the 'equatorial plane' (p = 7i"/2 are not flat. 

The virtue of these coordinates is that the form of the line element is invariant under 
the following transformation 

X C{x + X) 

y ^ C{y + Y) 

z ^ Cz (121) 

With a suitable choice of the constants X, Y and C we can place one of the two points to 
which the correlation function refers at an arbitrary position, while leaving unchanged the 
form of the line element and therefore the geodesic distance to the other point. 

If the point x',y',z' corresponds to the point r = in the spherical coordinate system, 
one can orientate the axes so that [|l^] 

X — x' = z' cos (/)sinhr/cosh(r/ — r) 
y — y' = 2;'sin(/)sinhr/cosh(r/ — r) 

z = z' coshr//cosh(T/ — r) (122) 

where tanhr^ = cos 9. 

Now consider the mode expansion. We look for eigenfunctions of the form 



^fcq^ = J dciiFk^Jz)e'^^-^ (123) 

where x is the vector with components x, y. Substituting this expression into the Laplacian 
gives a second order equation for -Ffcqj^. One of its two linearly independent solutions is 
1 17, 14 1 zKiq{q±z) where K is the modified Bessel function. This solution vanishes at 
z = oo, and up to normalisation is the only solution with that property. It also vanishes at 



z = (spatial infinity in opposite direction) for real k. Note that Kig is real for real q"^ and 
imaginary for imaginary q. 

We will define the sub-curvature modes by 

Wk^^ (x, z) = N{k)e'^^-^zKig{q^z) (124) 

where 

NHk) = i^S^ (125) 
The normalisation has been chosen to satisfy the orthonormality condition 

J T^fcq^(x, z)t^fc,q'Jx, ^)dV = <52(q^ - ql)5(g - q') (126) 

where the volume element is dV = dxdydz / z"^ . This condition is equivalent to the orthonor- 
mality of the functions (27r)~^e*'i^ '' plus the relation |73, 75| 

27r~2 poo 

sinh(7rg) / Kig{z)Kig.{z)dz/z = 5{q - q') (127) 



The expansion of a generic perturbation [^] in terms of these functions is 

/(x,z) = l^dqj dV/fcq^W^fcqJx,z) (128) 



Following Wilson [14|, we construct a Gaussian perturbation by assigning independent 
Gaussian distributions to the real and imaginary parts of the coefficients. Translation 
invariance in the x plane requires that the real and imaginary part of each coefficient has 
the same mean square, so we define the spectrum by (cf. Eq. (|2^ 



27r2 



The correlation function is 

roo r)—2 i\T2(h.\ I- 

C/(x,z,x',/)= J dq 2 + J d\^e^^^-^^-^hK,g{q^z)z'Kig{q^z') (130) 



It is invariant under Eq. (|121| ), because the factor C appearing in Eq. ( |121| ) can be absorbed 
into the definition of the integration variable q^ . 

Using Eq. (|122| ) and integrating over the angular direction in the qj^ plane, the correla- 
tion function becomes 



Jo 7r^[q^ + i) Jo 



(131) 



where 



From Eq. (8.13.30) of (7 



cosh(ry — r) 
sinh r 

1 = r 133 

cosh(77 — r) 



where 

7^ + + 1 

u = — = coshr (135) 

Ip 

This is a function only of r, and using Eq. (^) one finds that the correlation function is 
given by Eq. (p^). 

Now consider the super-curvature modes —1 < < 0. The normalisation factor N{k) 
becomes purely imaginary so it is convenient to drop the i factor, defining 

iV^(^ = J^i^^ (136) 

The super-curvature contributions are 

f'i^^z) = j\{iq)Jd\^fk^^Wk^Jx,z) (137) 
27r2 



The correlation function is given by Eq. (131) with q — > iq, and since Eq. (134) is valid for 
— 1 < ig < 1 this proves Eq. (|4^). 

The spherical and flat-surface expansions are equivalent, at least in the present context, 
because they both lead to a Gaussian perturbation with the same correlation function. 



The flat-space limit of the flat-surface expansion 

We end by looking at the flat-space limit of the flat-surface expansion. Though not strictly 
necessary for our purpose, it is extremely instructive, and does not seem to have been given 
before. 

The coordinates x, y and z become Cartesian in a a small region around x = y = and 
z = 1, and then qx, qy and q^ = ^Jq"^ — Qx ~ Qy the would-be components of the vector 
in the Fourier expansion. But the expansion Eq. ( |128|) does not restrict the range of qx and 
qy, so it will include both real and imaginary q^. In other words it will include hyperbolic 
functions as well as circular ones. Evaluating the limiting behaviour of Kiq confirms this. 
One finds ||7^ for real q^ 

Kigiq±z) ^ A sm{B + q,z) (139) 

where z = z — 1, A = \f^qz ^^^e"'^''/^ and i? = ^ -|- g cosh^"'^(g/g_L) — q^. For imaginary q^ 
one finds 

Kiq{qLz) ^ Cexp(Z) - \q,\z) (140) 

where C = (2|g^|/7r)i/2 and D = -q±sm-'^{q/\qz\). 

The second expression becomes infinite when \qz\z — > — oo, but even so a translation 
invariant correlation function will result when the two expressions are substituted into 
Eq. ( [1311 ). This example serves to remind us that we should take nothing for granted when 
considering which modes are allowable. 
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